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\0 ' Abstract 

o 

' We consider the equatorial circular motion of a test particle of specific 

£jl charge q/m <C 1 in the Kerr-Newman geometry of a rotating charged 

^ ^ 1 black hole. We find the particle's conserved energy and conserved 

bi)' projection of the angular momentum on the black hole's axis of rotation 

as corrections, in leading order of q/m, to the corresponding energy and 
angular momentum of a neutral particle. We determine the centripetal 
' force acting on the test particle and, consequently, we find a classical 

pseudo-Newtonian potential with which one can mimic this general 
relativistic problem. 
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In 1924, Manev (see also [2]) introduced a classical potential in order 
to modify the celestial mechanics in accordance with the general-relativistic 
description. To describe the motion of a particle of mass m in the static 
field of universal gravitation, due to mass M, Manev replaced the mass m 
with m = m,Q exp(M/r), where mo is an invariant and geometrized units 
G = 1 = c are used. This led to the following modification of Newton's 
gravitational law: 



where A and B are some real constants and M is the mass of the cen- 
tral body, was originally introduced by Newton himself in his Philosophiae 
Naturalis Principia Mathematica (Book I, Article IX, Proposition XLIV, 
Theorem XIV, Corollary II), in order to describe deviations in the Moon's 
orbit from the Keplerian laws. However, modification of Newtonian physics 
for mimicking general-relativistic results is attributed to Manev. We will 
refer to potential of the type (J2J) as to Manev potential. 
Modeling the particle motion in the most general stationary asymptotically- 
flat space-time — that of Kerr and Newman IU El — is a fundamental 
problem with many astrophysical applications. This area has been under in- 
tense investigation and various pseudo-Newtonian potentials for Kerr jHJ |5] 
and Schwarzschild [El IE] geometries have been proposed 0: interesting phys- 
ical applications, for example in satelite theory, have been studied in jS], 
different mathematical aspects in (Hj. In this paper we model the equatorial 
circular motion of a lightly charged particle {q/m <C 1) in Kerr-Newman 
geometry by introducing a classical pseudo-Newtonian potential. An imme- 
diate area of application of a classical pseudo-Newtonian effective potential 
of a charged particle is the description of accretion disks of charged particles 
around a rotating charged massive centre. In 1989, De Rujula, Glashow 
and Sarid jTU] proposed the idea that dark matter in our Universe is made 
of charged massive particles rather than neutral particles. These authors 
studied a galactic halo of such charged massive particles together with its 
cosmological implications, including galactogenesis. Gould et al. have 
shown that if charged massive particles (with masses bewteen 10 2 GeV and 
10 16 GeV) made up the dark halo of the Galaxy, then they would be present 
in large numbers in disk stars. Charged particles and their effects on galax- 
ies, stars and planets have been studied by Dimopoulos et al. |12j . 
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Potential of the type: 
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Rotating charged black holes are described by the Kerr-Newman geometry 
[21 [I] . The Kerr-Newman metric in Boyer-Lindquist coordinates ^Hl ( see 
also 5 ), in geometrized units G = 1 = c, is given by: 



ds l 
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(dt — a sin 2 6 
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In the above, M is the mass of the centre, a > — the specific angular 
momentum of the centre (i.e. angular momentum per unit mass) and Q - 
the charge of the centre. 

The motion of a particle of mass m and charge q in gravitational and elec- 
tromagnetic fields is governed by the Lagrangian: 



L 



1 

2 9lj 



dx l dx J 
dX dX 



q , dx l 
m dX 
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In the above, A is the proper time r per unit mass m: X = r/m and A is 
the vector electromagnetic potential, determined by the black hole's charge 
Q and specific angular momentum a: 



Ajdx 1 



Qr . o 

— 2~ \dt — a sin i 



P- 
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(The magnetic field is due to the dragging of the inertial reference frames 

into rotation around the black hole.) 

The equations of motion for the particle are: 



d?x l 



dxi dx 
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dr 2 3 dr dr m 3 dr ' 
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where F = dA is Maxwell's electromagnetic tensor and T l - k are the Christof- 
fel symbols. 

For Kerr-Newman geometry, the geodesic equations (jSJ) can be written as 
PH (see also [El): 



p 2 ^- = -a 2 E sin 2 6 + aJ + 
dX 



r 2 + a 2 
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where E = dL/dt is the conserved energy of the particle, J = dL/d<f> is the 
conserved projection of the particle's angular momentum on the axis of the 
black hole's rotation (dots denote differentiation with respect to A). K is 
another conserved quantity given by: 



K = V l + cos 2 9 a 2 (m 2 - E 2 ) + 



1 



sin 



(13) 



Here pg is the ^-component of the particle's four-momentum: pg = dL/d6. 
The radial equatorial (8 = vr/2) motion of the particle can be modelled as 
a one-dimensional motion of a classical particle in effective potential V{r) 
determined from (jl()|) : 



1 / dr \ 2 
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We now focus on circular orbits in the equatorial plane. For circular orbits, 
dr/dX = 0, both instantaneously and at all subsequent times (orbit at a 
perpetual turning point) [EllEj. This implies: 



V{r) 



0, 



(15) 



dV(r) 
dr 



0. 



(16) 



Equations ((T5|) and (fTB|) can be viewed as a system of simultaneous equations 
from which one can determine the conserved energy E of the particle and 
the conserved projection J of the particle's angular momentum on the axis 
of the black hole's rotation in terms of the parameters M, Q, a, q, and m. 
These equations are non-linear and we will not be looking for their general 
solution. Instead, we will assume that the particle's specific charge is small, 
i.e. q/m « 1, and we will expand E and J, as functions of q/m, with m 
fixed, in Taylor series near q/m = 0: 



E (—) = E(0) + E'(0)(^-) + 0[(-) 2 l = E + qE 1 + oU^-V] , (17) 
j(-)= J(0) + J'(0)(^)+O\(^) 2 } = J + qJ 1+ o\(^) 2 }. (18) 
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In the above expansions, one can identify E(0) = Eq as the conserved energy 
and J(0) = Jo as the conserved projection of the angular momentum of a 
neutral particle in Kerr-Newman Geometry. These are known ^ 



E r 2 - 2Mr + Q 2 ± ay/ Mr - Q 2 

m -< /- 2 - 3Mr + 2Q 2 ± 2ay/Mr - Q 2 



r\ r 



Jo _ ± a(Q 2 - 2Mr) + (a 2 + r 2 )yjMr - Q' 2 
'" ~ ,!.U/ • 2Q- : 2r/ v Mr Tp 



(19) 
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The upper sign corresponds to direct orbit (J > 0), while the lower sign 
corresponds to retrograde orbit (J < 0). 

The above expressions are valid for r > Q 2 /M. In addition, the condition 
M 2 > a 2 + Q 2 is necessary for the existence of a horizon (i.e. a black hole 
solution) in Kerr-Newman geometry . Orbits do not exist for all values of 
r > Q 2 /M. The question of existence and stability of bound circular orbits 
of neutral particles in Kerr-Newman geometry, as well as the thresholds 
for photon, marginally stable and marginally bound orbits are analysed in 
|18j (see also |16j for rotating black holes). Our goal is to find a classical 
pseudo-Newtonian potential for modelling the circular orbital motion of a 
lightly charged particle (q/m <C 1) in Kerr-Newman geometry by seeking 
expansions for large values of r (i.e. far from the black hole) which are above 
the threshold of a marginally stable orbit, i.e. within our approximation, all 
orbits are stable. 

Identifying the Keplerian angular momentum distribution J/E, the cen- 
tripetal force Fq acting on the neutral particle can be expressed as JJj : 
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Returning to the case of charged particle, we substitute the expansions (|T7|) 
and (JT8|) into equations ((T3|) and (fT^)) and we discard terms quadratic in 
q/m. We also use the fact that (fT9|) and (|20|) are solutions to equations (fT5|) 
and (|TB)) with q/m = 0. Thus, equations (|T5|) and ()16j) become: 

( r 2 + a 2 )El _ aJl _ a J °~l E ° (Jx - aEx) = Qr (22) 

(r z + a z )E - aJ 

2(M - r)J Ji - 2aM(J E 1 + E Ji) + 2(a 2 M + a 2 r + 2r 3 )E E 1 

= -aQJ + Q{a 2 + 3r 2 )E . (23) 
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These are now linear equations and they can be easily solved: 

£V(0) Q[± M Mr - Q 2 ) + ( r2 ~ ^Mr -a 2 + 3Q 2 )^/Mr - Q 2 



E, = 



m 



2r 
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As before, upper sign corresponds to direct orbit (J > 0), while lower sign 

corresponds to retrograde orbit (J < 0). 

The centripetal force acting on the charged particle is: 

2 
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The pseudo- Newtonian potential Vm(t) giving rise to this force is the Manev 
potential, 



V M (r) = - J Fdr, 
which, upon expansion over the powers of 1/r, results in: 
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From this expression it is evident that the leading term is reduced due to 
Coulomb repulsion if the two charges q and Q have the same sign. Opposite 
charges lead to additional attraction. By setting a = one obtains the 
effective pseudo-Newtonian potential for Reissner-Nordstr0m geometry |19j . 
and by further setting Q = — for Schwarzschild geometry (see also ) . 
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